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by Shih-Ying Wang^ and J. Marlin Smith 
Lewis Research Center 

SUMMARY 

Radial distributions of temperature were determined spectroscopically in a cesium- 
seeded, hydrogen-oxygen plasma jet. The plasma was generated at combustion cham- 
ber pressures ranging from 0. 5 to 2, 0 megapascals and for various seed ratios (1 to 
10 percent). The plasma was observed as the atmospheric exhaust from a NASA Lewis 
Mach 2 rocket test facility. 

Transverse profiles of the absolute integrated intensity were measured with the 
optically thin Csl lines (0. 5664 and 0. 5636 pm) at a range of axial positions downstream 
of the 5 -centimeter- diameter combustor -nozzle exit. Radial profiles of the emission 
coefficient were obtained from the measured transverse profiles of intensity by Abel 
inversion. Temperatures were then determined from the emission coefficients for con- 
ditions of local thermodynamic equilibrium using particle densities generated by a two- 
dimensional free jet computer program. Temperature results show the inherent ef- 
fects of compression and expansion pressure waves characteristic of a free jet exiting 
from a supersonic nozzle. 


INTRODUCTION 

Optical spectroscopy is frequently used to determine particle temperatures and 
densities in dense plasmas because it is one of the few diagnostic methods which does 
not disturb the plasma. The development of such optical diagnostics is therefore highly 
desirable for advanced combustion systems where physical probes cannot survive due 
to the severity of the environment, or where their sizes severely disturb the flow. 

This report discusses the spectroscopic measurement of emission intensities in a 
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cesium-seeded, hydrc^en- oxygen plasma free jet with a primary emphasis on deter- 
mining the radial temperature profiles in order to assess the spatial uniformity of in- 
jected seed. 


EXPERIMENTAL APPARATUS AND PROCEDURE 

A rocket test facility (ref. 1) at the NASA Lewis Research Center is being used to 
generate a combustion plasma. This plasma is a stoichiometric mixture of hydrogen 
and oxygen seeded with an aqueous solution of cesium hydroxide (CsOH) (75 percent 
CsOH, 25 percent H 2 O by weight) that is atomized in the oxygen flow line before enter- 
ing the combustion chamber. The plasma is accelerated in a Mach 2 nozzle and then 
expanded to atmosphere at the 5-centimeter-diameter exit. The sequence of gas injec- 
tion and ignition is controlled by a repeat cycle programmer to cover the running time 
from 5 to 20 seconds. Test- condition data such as flow rates and pressures are auto- 
matically sampled and processed by an IBM 360 computer; they are also recorded on a 
high speed oscillograph. Detection of the cesivun spectral lines is accomplished with a 
scanning monochromator having a reciprocal linear dispersion of 1. 6 nanometer per 
millimeter. Two optically thin Csl lines (0. 5664 and 0. 5636 pm) and a continuum at 
0. 5704 micrometer (for continuum correction) are found appropriate for absolute inten- 
sity measurement (see appendix A). 

The optical sjrstem is shown in figure 1. Plasma radiation first strikes the sta- 
tionary mirror, which is inclined at a 45° angle to both the optical and vertical axes. 

A second mirror is mounted below the stationary mirror on a motor-driven rotating 
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Spectrograph CS-78782 


Figure 1. - Orientation of spectrograph ic equipment. 
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base in such a fashion as to be also inclined at a 45° angle to both the optical and ver- 
tical axes. Thus, in traversing the mirror system, the radiation undergoes two 90° 
rotations, one about the vertical axis and one about the optical axes. Thus, the plasma 
image is focused on the spectrometer entrance slit such that the plasma x-axis is 
parallel to the slit. With reference to figure 2, the slab of plasma with radiant inten- 
sity Qy therefore has a width Ay defined by the spectrometer slit width. A fishtail 
diaphragm at the entrance slit is used to define the plasma length Ax irradiating the 
spectrometer optics. 

The rotating-mirror platform scans the plasma image across the slit and thereby 
provides a transverse profile Qy(y). The photomultiplier response to this transverse 
scan is recorded on both a memory oscilloscope and the high speed visicorder. System 




(b) Transverse profile of intensity. 
Figure 1 . - Plasma profile. 
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calibration for absolute intensity measurements is accomplished with a tungsten ribbon 
filament lamp mounted in place of the plasma on the plasma axi^. 

Spectrometer entrance slit widths are characteristically 30 to 60 micrometers, 
with corresponding exit slit widths being 600 to 1200 micrometers. 


SPECTROSCOPIC METHOD 

A gas is said to be in a condition of local thermodynamic equilibrium (LTE) when 
the electrons, atoms, molecules, and ions are in equilibrium with themselves, but not 
with the photons. In order that the plasma be in LTE, it is necessary that the collisional 
rates which populate and depopulate the various energy levels of the plasma species ex- 
ceed the corresponding radiative rates. Cool (ref. 2) has shown that LTE exists in an 
atmospheric potassium -seeded plasma for an electron temperature of 2500 K and an 
electron number density of 10 per cubic centimeter. In this case the rates of radia- 
tive deexcitation and recombination were well below that of electron collisional deexcita- 
tion and recombination. These conditions approximately represent the conditions of the 
present experiment, except that the major plasma constituent is water vapor rather than 
air. However, the radiative/collisional rate ratio for water is less than for air, thereby 
further justifying the assumption that the present plasma is in LTE. 

Line emission is determined by bound-bound transitions between upper level m and 
lower level n. The radiated power per unit voliune, per unit solid angle, and per unit 
wavelength interval is (ref. 3) 
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where L^ is taken to be the normalized Lorentz profile from reference 4: 
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for a line center at Xq with halfwidth ■^nm transition probability between 

lower level n and upper level m, and is the number of density at upper level m 

related to the total number density N of neutral cesium atoms by the Boltzmann re- 
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Relating the transition probability and the absorption oscillator strength 

by 
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gives the emission coefficient radiated over the bandwidth AA. of a spectral line as 
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where 



(6) 


The function F is plotted in reference 5 as a function of the volumetric absorption 
coefficient pKu and the ratio of bandwidth to halfwidth AA/Wj^y 2 approximately 

equal to one for the present experiment. The number density of neutral cesium atoms 
N„„, obtained from the ideal gas law, is given by 

oil 

= 7. 24X10^^ £§ (7) 

sn ^ 

where p is the gas pressure in pascals, T the temperature in kelvins, and S the mole 
fraction of neutral cesium atoms as calculated by a chemical equilibrium composition-' 
program (ref. 6). Values of the energy of the upper state E^, the statistical weight at 
the lower state g^^, and the oscillator strength are found in references 7 and 8. 
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ABSOLUTE INTENSITY CALIBRATION 


As mentioned previously, the optical system is caUbrated using a tungsten ribbon 
filament lamp moimted such that the filament is on the plasma axis. With the spectom- 
eter set at the wavelength of the cesium lines, a series of measurements are made at 
increasing filament brightness; each measurement consists of a photomultiplier signal 
voltage and a f ilam ent brightness temperature determined with an optical pyrometer. 

It is desired to formulate a spectrometric calibration function which will relate the 
photomultiplier signal voltage to Q , the radiant power flux per unit solid angle 

within the bandwidth AX of the spectrometer at the test wavelength X^. Recognizing 
that for calibration purposes the test radiation is rising from the tungsten filament, Q 

s 

is given by 


= r(xpe(Xj, T^Qp(Xt, T^)AX (8) 

where t(X^ is the transmittance of the lamp glass envelope, e(X^, T^ the spectral emis- 
sivity of the tungsten filament, Qp(^^, T^) the Wien form of the blackbody radiation func- 
tion, and AX the bandwidth of the spectrometer. 

The true temperature T^, of the tungsten filament is obtained from its brightness 
temperature observed by a pyrometer of effective wavelength Xp. Since the py- 
rometer is calibrated in terms of a blackbody, it follows that 


Qp(Xp,T^) = r(Xp)€(Xp.T,)Qp(Xp,T,) 


from which it follows that 



— + — ln[r(Oe(X ,T )] 

rp Q y P 


(9) 


( 10 ) 


Thus, from equations (8) and (10) values of Q are obtained for each brightness temper- 
ature measurement which, together with the corresponding Vp^ values, yields the 
desired spectrometric calibration function by a least squares fit in the form 

Q 5 (X() . Co t + OgV^m (11) 
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ABEL INVERSION 


Generally 18 to 32 data points (Qy, y) are obtained from a transverse scan of a 
plasma profile (fig. 2(a)). Following the method of Cremers and Birkebak (ref. 9), 
these data are divided into three zones and a least square fourth-degree polynomial is 
fitted to each zone. The zones overlap one another by four data points, and the criterion 
for selection of the intersection between the zones is based on which intersection yields 
the best fit to the data. 

Anal 3 dically, the Abel transformation of to radial emission coefficients is 
accomplished by the relation 



For digital computer processing, however, the numerical technique of Nester and Olsen 
(ref. 10) is employed. In this technique both the radius and the y-axis are divided into 
increments "a” such that 


rj^ = ka (13) 

yj^ = na (14) 


and 


v = Na 
•'max ■ 


(15) 


where y is the value of y where = 0 as shown in figure 2(b). From the geom- 

1X13JC 

etry of the inversion process shown in figure 2(a), it is seen that y„_^ is numerically 
equal to R. 

With these identifications, equation (12) may be manipulated into the form 
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where 


Qlj = Qj. r' = ka 
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TEMPERATURE PROFILE RATIOCINATION 

From equations (5) and (7) it can be seen that the emission coefficient is functionally 
dependent on temperature, plasma static pressure, and the mole fraction of neutral 
cesium atoms: 


Q(r)~EiEM£2exp[-^l 

T(r) LkT(r)J 

Obviously, the radial temperature profile can be determined from Q(r) only if p(r) 
and S(r) are known from other sources. For this work, p(r) was obtained from a two- 
dimensional supersonic free jet computer program which uses the method of character- 
istics (ref. 11 and private communication with A. R. Bishop, NASA Lewis Research 
Center) to relate pressure to temperature, input composition, and nozzle geometry. 
S(r) was obtained from a chemical equilibrium program (ref. 6) which thermodynami- 
cally relates S(r) to T(r) and input composition. The relations p^.(T^, composition, 
geometry), S^(Tj,, composition), and Q^(T^, p^, sp were iteratively solved for that 
T^ yielding a self-consistent solution. 


RESULTS AND DISCUSSION 

A set of profiles in sequence of integrated intensity, emission coefficient, pressure, 
and temperature is shown iii figures 3 to 6, respectively, for combustion chamber pres- 


8 



Emission coefficient. W/m^-sr Intensity. W/m^-sr 



y-Axis, cm 


Figure 3. - Integrated intensity as function 
of y-axis at combustion pressure of 10^ 
newtons per square meter. 



Figure 4. - Emission coefficient as function 
of radius at combustion pressure of 10” 
newtons per square meter. 




Figure 5. - Pressure as function of radius at 
combustion pressure of 10“ newtons per 
square meter. 



Figure 6. - Temperature as function of radius 
at combustion pressure of 10“ newtons per 
square meter. 
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sure of 1. 0 megapascal. Each figure contains five curves representing data taken at 
different axial positions downstream of the nozzle exit. The plume nodes occur at po- 
sitions 3. 8 and 19. 0 centimeters, approximately. The shapes of the temperature pro- 
files in figure 6 reveal the inherent effects of compression and expansion pressure waves 
characteristic of a supersonic free jet. The centerline temperatures vary from about 
2900 to 3200 K. The free boundary contour of the plasma pliune has 2. 4 to 2. 8 centi- 
meter radii as defined by a temperature of about 2400 K, below this the intensity is too 
small (<1. 0 W/m - sr) to record. 

The effect on the temperature profile of lowering the chamber pressure from 1. 0 to 
0. 5 megapascal is shown in figure 7. The profile is generally steeper and the plume 
size is smaller (radii, 1. 8 to 2. 1 cm). The corresponding pressure profile is shown in 



Figure 7. - Temperature as (uncljon of radius 
at combustion pressure of 5x10^ newtons per 
square meter. 
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figure 8. Raising the chamber pressure to 2. 0 megapascals results in the plume radii 
increasing to 2. 7 to 3. 7 centimeters, as shown in figure 9, while some of the tempera- 
ture profiles decrease from their maximum of about 3000 K to approximately 2600 K at 
the centerline. The corresponding pressure profile is shown in figure 10. 

The effect of seed fraction can be seen from figure 11 for the same chamber pres- 
sure and axial position. The rate of change of peak intensity decreases with increasing 
seed fraction. 

Temperature profiles were also measured for some of the previous operating con- 
ditions using the Csl 0. 5636 -micrometer transition. One such typical measurement is 
shown in figure 12. On the average, temperature measurement using the 0. 5636- 
micrometer line agree to within 4 percent with those obtained using the 0. 5664- 
micrometer line. 

The experimental uncertainty in temperature due to intensity measurement, data 
reduction, and calibration procedures is estimated to be about 5 percent. 


2.0x105 



Radius, cm 

Figure 8. - Pressure as functton of radius at 
combustion pressure of SxlCr newtons per 
square meter. 
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Radius, cm 

Figure 9. - Temperature as function of radius at combustion 
pressure of 2(10" newtons per square meter. 



0 12 3 4 

Radius, cm 


Figure 10. - Pressure as function of radius at combustion 
pressure of 2(10" newtons per square meter. 
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0 .02 .04 .06 . 08 .10 


Seed fraction 

Figure 11. - Peak intensity ^ function of seed fraction at 
combustion pressure of 10° newtons per square meter 
and X = 5.4 centimeters. 


Run Seed Wavelength, 

fraction, X, 

3100| — percent pm 


382 5.2 0.5665 

372 5.0 . 5636 



Radius, cm 


Figure 12 - Radial temperature profiles at 
0. 5664 and 0. 5636 micrometer wave- 
lengths, combustion pressure of 10° 
newtons per square meter, and x = 

3.8 centimeters. 
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CONCLUDING REMARKS 


Results of the spatially resolved temperatures have demonstrated the feasibility of 
the diagnostic technique of absolute intensity measurement for the seeded combustion 
plasma jet, if the density profile is known. The work provides a tool to study the ther- 
mal field of an axisymmetric MHD duct. 

If the plasma were measured from side and top (i. e. , at different angles in the 
plane perpendicular to the flow axis), the axisymmetry would be better assessed. 

The reader may ask why temperature was not measimed using a two lines method 
so that density might be determined using T(r) as the independently measured quantity. 
This was tried, but the fluctuations in the temperatures so obtained were intolerably 
large. These fluctuations may have arisen from two sources: (1) the errors in the 
Abel inversion technique, or (2) the fact that sequential rather than simultaneous scans 
were measured due to instrumental limitations. The failure of this initial approach 
dictated the development of the semi -theoretical approach herein described. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 
506-25. 
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APPENDIX A 


METHOD FOR ESTIMATION OF THE OPTICAL THICKNESS 


The optical thickness estimation has been reported before (refs. 12 and 13) by cal- 
culating the volumetric absorption coefficient: 


pKv 


e^ ^mn ^n 


N. 


sn 
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where the line halfwidth is dominated by Lorentz broadening. This pressure broadenii^ 
by different gas particles is given in references 14 and 15 by 


fN 

c ^opt 

7T ^ 



(A2) 


An experimental value for the optical cross section Q , of the Csl 0, 8521- 

^ 1/5 

micrometer line is given in reference 15 in the form Qopt/"^g^ • optical cross 
sections of other Csl lines can be obtained from this using the following relation from 
reference 16: 


, Qopt 

o 

0.8521 pm 0.8521 pm 



(A3) 


where C is the van der Waals constant of interaction. The quantum mechanical average 
of the square of the so-called "radius” of the initial (a = i) and final {a = f) energy 
states is given in reference 3 as 




+ 1 - 3i (I 
a a 



(A4) 


where Ett = 13. 595 electron volts, E = 3. 893 electron volts, and I is the orbital 

Xl f ao ■’ 01. 

quantum number of the state i or f. 

For a gas temperature of 3000 K, the cross sections of the Csl, 0. 5664 and 0. 5636 
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micrometer lines are 6. 2 and 7. 9 square nanometers, respectively. For = 0. 001, 
pKu < 0. 15 per centimeter for both lines, which are considered approximately optically 
thin. 
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APPENDIX B 


SYMBOLS 


■^nih 

^0 

®k, n 
C 

c 


V *^2’ ^3 
^4 

k’ n’%n-l 
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e 
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inn 

^m’ ^n 

go 
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M_ 


M 


sn 


m. 
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N, 

N 
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transition probability 

constant (ref. 3, p. 99), first Bohr radius, 0. 52917 A 

defined in ref. 9 

Van der Waals constant 

speed of light 

defined by eq. (11) 

1.44x 10"2 (m)-(K) 
defined by eqs. (14) to (16) 

first ionization potentials of H and Cs atoms, respectively 
energy of upper level 
electron charge 

ratio of spectral line energy in bandwidth AX to total spectral 
line energy 

absorption oscillation strength 

statistical weight of upper and lower state 

statistical weight of ground state 

Planck’ s constant 

Boltzman’s constant 

Lorentz broadening spectral line profile 

orbital quantum number of state a 

molecular weight of gas 

atomic weight of Cs seed 

electron mass 

number density of gas particles 
number density of Cs atoms at upper level 
number density of neutral Cs atoms 
pressure 
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^opt 

% 

Qr 

Qs 

Qy 


2 2 


Sf 

T 




X 


y 

^0 

X 

p 

AX 

pKu 

T 


optical cross section 
Planck function 
emission coefficient 
defined by eq, (8) 
integrated intensity 
radius 

universal gas constant 
classical radius of electron 

square of classical electron radius of states i and 

mole fraction of neutral Cs atoms 

temperature 

brightness temperature 

true temperature 

photomultiplier tube voltage output 

spectral line width at half height 

axial position downstream of nozzle exit 

axis perpendicular to x 

partition function of Cs atom (=2) 

emission coefficient 

wavelength at line center 

effective wavelength of pyrometer 

bandwidth 

volumetric absorption coefficient 
transmittance 

collisional broadening halfwidth in wave numbers 



REFERENCES 


1. Smith, J. M. ; Results and Progress on the NASA Lewis H 2 -O 2 MHD Prc^ram. 

15th Symposium on the Engineerii^ Aspects of Magnetohydrodynamics, Pennsyl- 
vania Univ. , 1976, pp. DC. 2-DC. 6. 

2. Cool, Terrill Alan; Recombination, Ionization, and Nonequilibrium Electrical Con- 

ductivity in Seeded Plasmas. Ph. D. Thesis, California Institute of Technology, 
1965. 

3. Griem, Hans R. ; Plasma Spectroscopy. McGraw-Hill Book Co. , Inc. , 1964, p. 99. 

4. Mitchner, M. , and Kruger, Charles H. , Jr. : Partially Ionized Gases. John 

Wiley & Sons, Inc. , 1973 . 

5. Morris, J. C . ; Krey, R. U . ; and Garrison, R. L. : Radiation Studies of Arc 

Heated Nitrogen, Oxygen, and Argon Plasmas. Tech. Rep., Ijunel965- 
15 June 1967, Avco Corp. , 1968. (AFARL-68-0103, AD-676651.) 

6. Gordon, Sanford; and McBride, Bonnie J. : Computer Prc^ram for Calculations of 

Complex Chemical Equilibrium Compositions, Rocket Performance, Incident and 
Reflected Shocks, and Chapman-Jouquet Detonations. NASA SP-273, 1967. 

7. Agnew, L. ; and Summers, C. : Quantitative Spectroscopy of Cesium Plasmas. 

Proceedings of the 7 International Conference on Phenomena in Ionized Gases, 
Vol. n, B. Perovic and D. Tosic, eds. , Boris Kidric Inst. Nuclear Sciences 
(Giograd), 1966, pp. 574-580. 

8. Moore, Charlotte E. : Atomic Energy Levels as Derived from the Analyses of 

Optical Spectra. NSRDS-NBS 35, Vol. 3, National Bureau Standards, 1971. 

9. Cremers, Clifford J. ; and Birkebak, Richard C. : Application of the Abel Integral 

Equation to Spectrographic Data. Appl. Opt., vol. 5, no. 6, June 1966, pp. 1057- 
1064. 

10. Nester, O. H. ; and Olsen, H. N. r Numerical Methods for Reducing Line and Sur- 

face Probe Data. SIAM Rev. , vol. 2, no. 3, July 1960, pp. 200-207. 

11. Hayes, Wallace D. ; and Probstein, Ronald F. : Applied Mathematics and Mechan- 

ics. . Vol. 5, Hypersonic Flow Theory, Academic Press, 1959, chap. 7. 1. 

12. Wang, S. Y. ; and Benenson, D. M. : Radial Distribution of Electron and Gas Tem- 

peratures in a Potassium-Seeded Argon Plasma. International Conference on 
Plasma Science, IEEE, 1976, Session 2B11, p. 49. 

13. Wang, Shih-Ying: Experimental Investigation of a Potassium -Seeded Argon Plasma. 

Ph.D. Thesis, State Univ. of New York at Buffalo, 1975. 


20 



14. Mitchell, Allan C. G. ; and Zemansky, Mark W. : Resonance Radiation and Excited 

Atoms, Univ. Press, Cambridge, 1961, p. 170. 

15. Hiram oto, Tatsumi: Rates of Total and Local Radiative Energy Losses in Non- 

equilibrium Plasmas. J. Phys. Soc. Jpn. , vol. 26, no. 3, Mar. 1969, pp, 785- 
801. 

16. McGregor, Douglas David: Electronic Nonequilibrium in a Supersonic Expansion of 

Ionized Gas. Appendix M, Optical Cross Sections for the 5P-3D Transition. 
SU-IPR-474, Stanford Univ. , 1962, pp. 247-248. (AFOSR-72-2049TR, AD- 
750943. ) 


21 


1. Report No. 2. Government Accession No. 

NASA TP 1162 

3. Recipient's Catalog No. 

4. Title and Subtitle 

TEMPERATURE DISTRIBUTIONS OF A CESIUM-SEEDED 
HYDROGEN -OXYGEN SUPERSONIC FREE JET 

5. Report Date 

February 1978 

6. Performing Organization Code 

7. Author(s) 

Shih-Ying Wang and J. Marlin Smith 

8. Performing Organization Report No. 

E-9267 

10. Work Unit No. 

506-25 

9. Performing Organization Name and Address 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Paper 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

Radial distributions of temperature were determined spectroscopically in a Cs-seeded, H 2 -O 2 
plasma jet. The plasma was generated at combustion chamber pressures ranging from 0. 5 to 
2. 0 megapascals and for various seed ratios (1 to 10 percent). The plasma was observed as 
the atmospheric exhaust from a NASA Lewis Mach 2 rocket test facility. Transverse profiles 
of the absolute integrated intensity were measured with the optically thin Csl lines (0. 5664 and 
0. 5636 pm) at a range of axial positions downstream of the 5-cm-diam combustor-nozzle exit. 
Radial profiles of the emission coefficient were obtained from the measured transverse pro- 
files of intensity by Abel inversion. Temperatures were then determined from the emission 
coefficients for conditions of local thermodynamic equilibrium using particle densities gen- 
erated by a two-dimensional free jet computer program. Temperature results show the in- 
herent effects of compression and expansion pressure waves characteristic of a free jet exit- 


ing from a supersonic nozzle. 


17. Key Words (Suggested by Author(s) ) 

Plasma physics 
Spectroscopy 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 75 

19. Security Classif. (of this report) 

Unclassified' 

20. Security Classif. (of this page) 
Unclassified 

21. No. of Pages 
22 

1 

22. Price* 


For sale by (he National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1978 





















National Aeronautics and 
Space Administration* 

Washington, D.C. 

.20546 

Official Business: 

Penalty for Private Use, $300 


SPECIAL FOURTH CLASS MAIL 
BOOK , 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 V 



, 


rUASA 


N. 




U ■ 


POSTMASTER 


■ / . 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 




^ .. 






. c 




J ^ \ ■ , 

' At* // 'p 




. . . 




V 




/ 


. 7 

< N 


A 


■ A- ' \ 




■ t 




